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ABSTRACT 

A finite element analysis of the Iosipescu shear test for unidirectional and cross-ply 
composites is presented. It is shown that an iterative analysis procedure must be used to 
model the fixture-specimen kinematics. The correction factors which are needed to 
compensate for the nonuniformity of stress distribution in calculating shear modulus are 
shown to be dependent on the material orthotropic ratio and the finite element loading 
models. Test section strain distributions representative of typical graphite-epoxy 
specimen are also presented. 

Keywords: Numerical Analysis, Composite Materials, Iosipescu Specimen, Shear 
Modulus, Correction Factors 


INTRODUCTION 


The development and evaluation of in-plane shear test methods for the 
determination of shear modulus and strength has received considerable experimental and 
analytical attention from the composite materials industry, [1-5]. An ideal shear test 
method would produce a pure and uniform shear stress field in the test section of the 
specimen throughout the linear and nonlinear response regimes. Although the Iosipescu 
shear test method does not possess these ideal characteristics, it has gained wide 
acceptance for measuring the in-plane shear response of composite materials. 
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Two versions of the Iosipescu shear fixture and specimen configuration [4, 6] have 
been developed for composite materials, as depicted in Figs. 1 and 2. The original 
fixture (Wl) produced a very small region of uniform shear stress and significant normal 
strains in the specimen test section. To overcome these deficiencies a modified fixture 
(W2) and specimen configuration [6] were developed. The modified fixture and 
specimen produced a larger test section, larger fixture-to-specimen contact regions and 
the innermost fixture-to-specimen contact points were moved farther away from the test 
section. The uniform shear stress region increased but normal strains still exist in the test 

section [7]. 

A considerable amount of experimentation and numerical analysis of the Iosipescu 
shear test have been performed [6-19]. Numerical studies have proven useful in 
understanding the influence on shear response of notch angle and radius [6], ply 
orientation [6,7,1 1,13,16] as well as failure mechanisms and failure location [13], 
uniformity of stress-strain distribution in the test section [7-9 ,11-16] and the 
development of shear modulus correction factors [7,1 1,16]. In any numerical study it is 
essential to represent accurately the kinematics of the fixture and specimen because 
boundary conditions and load introduction methods can significantly influence the stress 
state in the test section and the predicted failure location. 

Unfortunately, there seems to be a lack of consensus of how to represent the 
boundary conditions and load introduction in the Iosipescu specimen. For example, 
Barnes et al. [13] used uniformly distributed applied forces, Sullivan et al. [9] used 
concentrated forces and Walrath and Adams [8] used prescribed vertical (y) 
displacements at the specimen-fixture contact regions in their finite element analyses of 
the Wl test specimen. The use of concentrated forces can cause significant local 
deformations in the specimen and are unrepresentative of the experimentally observed 
deformations. A uniform applied force assumes the applied force distribution is known a 
priori. In recent finite element analyses of the W2 specimen by Adams and Walrath [6] 
and Pindera et al. [11,16], prescribed vertical displacement loading conditions were used 
but different types of constraints (hinges, u=v=0, and rollers, v=0) were applied to the 
stationary part of the specimen-fixture contact regions. Furthermore, the same loading 
condition was applied to the 0° and 90° specimens [11,16]. 

The objective of this study is to conduct a linear-elastic finite element analysis of 
the modified Iosipescu shear specimen using a more realistic method of modeling the 
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load transfer between the fixture and specimen and the displacement conditions at the 
specimen/fixture interface. The analysis is based upon an iterative scheme without any 
bias on the load distribution into the specimen. The reacdon force distribution along the 
fixture-to-specimen contact region, the strain state in the specimen and the shear modulus 
correction factors are evaluated. Three fiber orientations ( 0°, 90° and 0790°) are used in 
the analysis along with three materials having different orthotropic ratios (typical values 
of graphite-, Kevlar- and glass-epoxy). 


FINITE ELEMENT MODELS 


In the modified Iosipescu shear test, the left end section of the specimen is clamped 
into the stationary part of the fixture and the right end section of the specimen is clamped 
into the movable part of the fixture, see Fig. 2a. Assuming that the clamping process 
does not introduce any significant strains in the specimen, load is transferred only from 
the fixture to the specimen along the upper fixture-to-specimen contact region of the 
movable portion of the fixture. The lower fixture-to-specimen contact region of the 
movable portion of the fixture cannot "pull" on the specimen because it is not 
mechanically attached to the specimen. That is, tensile forces cannot be applied along the 
lower edge of the specimen by the fixture. However, a reaction force is developed 
between the specimen and fixture along the lower edge of the movable part of the fixture 
because the lower contact region of the fixture prohibits rotation of the specimen keeping 
the upper and lower surfaces parallel and horizontal, as depicted in Fig. 3. 

Local in-plane bending of the specimen occurs as load is applied to the specimen, as 
depicted in Fig. 3b. Along the outboard region of the upper fixture-to-specimen contact 
region and along the inboard region of the lower fixture-to-specimen contact region 
(movable portion of the fixture), the specimen can progressively separate from the fixture 
as a function of the applied load, see Fig. 3b, reducing the length of the contact region. A 
similar set of deformations occurs on the fixed portion of the specimen. 

In this study, the distance from the notch axis to the inner most load point was 
7.0mm. The mechanical properties of the three materials are presented in Table 1. The 
finite element model had 1468 membrane finite elements, see Fig. 4. Close to the notch 
root, the elements are refined such that the reduction of strains to near zero at the free 
boundary are achieved. Prescribed displacements, v=8= -0.05mm, were applied to the 
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fixture-to-specimen contact regions on the movable part of the fixture. The ABAQUS 
and PA TRAN codes were used for calculations and pre-, post processing, respectively. 

The sequence of steps for application of boundary conditions and applied 
displacements in this study are as follows: 

1 Displacement constraints (v=0 or u=v=0) to eliminate rigid body motion of the 
specimen are applied along the fixture-to-specimen contact region on the stationary (left) 
part of the specimen, see Fig. 5. 

2. Uniformly applied displacements are applied along the fixture-to-specimen contact 
region on the movable (right) part of the specimen, see Fig. 5. 

3 Reaction forces are calculated at the nodes where the applied and constrained 
displacements are introduced. Displacements along the fixture-to-specimen contact 
region are also calculated. 

4. If any reaction force were tensile or if the specimen interfered with the fixture 
(displacement interference) then the displacement constraint or the applied displacement 
is changed or removed from that finite element node. 

5. The analysis is performed with an updated set of applied displacements and 
boundary condition and the process is repeated starting with step 3 until convergence 
occurs. 

6. When convergence occurs the load applied to the specimen is the sum of the 
reaction forces on the movable (right) portion of the specimen. 


RESULTS AND DISCUSSION 

Force Distribution Alonn The Snecimen-to-fix tiire Contact Region 

Solution convergence usually occurred within 3 or 4 iterations. The resulting 
displacement and reaction force distribution was different for each of the 0°, 90° and 
0°/90° specimens, as depicted in Fig. 6. These force distributions are a function of the 
ply orientation and are representative of neither the concentrated nor uniformly 
distributed force conditions applied by others [9-13]. The load distribution along the 
upper fixture-to-specimen contact region on the movable part of the fixture is highly 
nonlinear with the highest forces closest to the center of the specimen. The force 
distribution along the upper surface of the movable part is not the same as the force 
distribution along the lower surface. The force distribution at the innermost fixture-to- 
specimen contact region on the 90° specimen is narrower than for the other ply 
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orientations because in-plane bending occurs causing the specimen to pull away from the 
fixture. The force distribution of the 0790° specimen is approximately the average of the 
0° and 90° specimens. 


Strain State In The Specimen 


Test section 


The normal and shear strains, normalized with respect to the average shear strains 
across the test section, are shown in Figs. 7-9 for 0°, 90° and 0790° specimens, 
respectively. The horizontal normal strains, £ x , are negligible for all 0°, 90° and 0790° 
specimens because no horizontal constraints are applied which is consistent with 
experimentally observed full-field results [7]. The compressive transverse normal strains, 
e are negligible for 90° and 0790° specimens but are significant for the 0° specimens. 

As shown in Fig. 7 for the 0° specimens, e y reaches a maximum at the center of the 
specimen and decreases to zero at the notch root. The £ y strain distribution between 
notch tips (at x=0) for the 0° specimen is approximately a parabolic shape. At the center 
of the specimen, £ y is 24.4 percent of y avg for graphite-epoxy and is 9.5 percent and 12.3 
percent of y av g for the Kevlar-epoxy and glass-epoxy, respectively. The compressive 
forces at the inner most fixture-to-specimen contact region creates £ y strains in the test 
section of the 0° specimens, see Fig. 7. The £ y strains are created in the test section 
because the 0° specimens bend and compress due to the low transverse stiffness of the 
laminate. The transverse compression of the 0° specimens decreases with increasing 
distance from the inner most specimen-to-fixture contact region. The 90° and 0790° 
specimens experience negligible £ y strains, see Figs. 8 and 9, in the test section because 

the 90° fibers resists the transverse compression of the specimen. 

The shear strain contours of the 0°, 90° and 0790° graphite-epoxy specimens are 
shown in Fig. 10 whereas the normal strain contours are depicted in Figs. 11 and 12. In 
the test section of the specimens where the strain gages would be attached, the shear and 
normal strains are approximately uniform. Strain contours for 0790° specimen have 
characteristics of both 0° and 90° specimens and its response is essentially the average of 
the response of the 0° and 90° specimens. 
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Fixture-to-specimen contact region 


In the 0° specimen, £ y is largest at the inner most contact region and decreases 
radially from the innermost contact region into the test section, Fig. 1 la. The E y in the 
test section is about 16% of the £ y at the innermost contact region. As shown in Fig. 12a, 
e x is small, about 14% of the £ y , at the innermost contact region. Shear strain y xy is 
similar in magnitude to £ y . For the 90° specimen, the largest £ y occurs at the two 
antisymmetric inner contact regions next to the notch root and only decreases vertically. 
Fig. lib. In the test section, there is no presence of £ y . The magnitude of £ x at the 
innermost contact region is compressive and is of the same magnitude as £ y and y xy . For 
0790° specimens, the longitudinal strain £ x contour is more representative of the strain 
contour of the 0° specimen whereas the £ y contour is similar to that of 90° specimen. 


Notch region 

Ideally, the specimen should fail in shear at the minimum cross section so that the 
notches should not provide any stress raising effect. Experimental observations [7,15] 
indicate that the ideal condition is not achieved and that the failure locations for the 0° 
and 90° specimens are found at regions close to the intersection of the notch root and 
notch flank, see Fig. 2b, rather than at the minimum cross section. The linear elastic 
finite element results obtained here indicate that a complex strain state occurs in the 
region where initial failure is generally observed. At the intersection of the notch root 
and flank the transverse normal strain, £ y , is tensile with a value of approximately 24%, 
see Figs. 10a and 1 la, of the largest shear strain y xy in the test section, and the shear 
strain y xy at the notch root-flank intersection is approximately 74% of the largest shear 
strain in the test section, for the 0° graphite-epoxy specimen. 

In the case of the 90° specimen, Fig. 12b, the longitudinal strain £ x is tensile, at the 

notch root- flank intersections further from the applied loading (point x in Fig. 12b), with 
a value of approximately 1 8% of the shear strain y xy at the center of the test section. The 
corresponding shear strain at the notch root-flank intersection is about 51% of the y xy at 
the center of the test section. 

It is beyond the scope of this linear elastic analysis to estimate the failure mode or 
strength of the Iosipescu specimens because of the significant nonlinear material and 
geometric effects that occur. However, this analysis will provide insight into failure 
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location and some of the important material and geometrical parameters that control the 
response of the specimen. The results from the finite element analysis indicate mixed 
mode failures are feasible in the notch root-flank region. This is supported by 
experiments [7,15] which indicate that the notch root-flank intersections are critical 
regions where failures are frequently initiated. Therefore, the use of the 0° and 90° 
Iosipescu specimens to measure shear strengths is ill-advised because the strength is not 
solely a function of the shear properties of the material. An additional complexity occurs 
for the 90° and 0790° specimens, that is the specimens twist [7] which produces mixed 
mode failures. 


Correction Factor 

The shear strain between notch tips is not constant for the different ply orientations 
and materials evaluated, see Figs. 13a- 13c. If a shear modulus G 12 were calculated using 
the experimentally measured shear strain (using strain gages), erroneous results would be 
obtained. To account for the nonuniformity of shear strain distribution in shear modulus 
calculation, correction factors are developed. Correction factors are obtained by dividing 
the shear strains, Ygage^ the center of the test section of the specimen by the average 
shear strain, Yavg. between the notch tips. The average shear strain, Ygage. can be 
analytically determined by averaging the shear strain over an area that would be covered 
by a strain gage. The correction factor and correct shear modulus are calculated by 

A 

7avg = h J k Txv 
2 


CF = SfBl 

Yavg 

q _ T avg 

^ Yavg Ygage (3) 

where h is the distance between two notch tips and G* is the apparent shear modulus 
(^avg/Ygage ) _ a comparison of shear strain distributions along the notch axis normalized 
with respect to the average shear strains for 0°, 90° and 0790° graphite-epoxy, Kevlar- 
epoxy and glass-epoxy specimens is shown in Fig. 13. The corresponding correction 


dy 

(1) 

( 2 ) 

Tav s ~ ri? _ rr v n * 
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factors are shown in Fig. 13 and depend on the material orthotropic ratios. For the 0° 
specimen, the correction factor is less than 1.0. The larger the material orthotropic ratio, 
the smaller the correction factor. For 90° specimens, the reverse trends were observed. 
Approximately, the correction factors can be expressed as, 

CF = 1.036 - 0.125 x log (§*■) 

(4) 

where E x and Ey are extensional stiffnesses in the longitudinal and transverse directions, 
respectively. Figure 14 is a plot of the correction factors as a function of material 
orthotropic ratio. Note that for the 0790° specimens, which have unit orthotropic ratio, 
the shear strain distributions are more uniform in the gage section. In the calculation of 
shear modulus, however, a correction factor is still necessary. The correction factor for 
the 0790° specimens is approximately the average of the correction factors of the 0 and 
90° specimens. Note that the correction factors depend on the material orthotropic ratios 
and are insensitive to the shear moduli to be measured. In Fig. 15, it is shown that 
increasing or decreasing the shear modulus of the 0° graphite-epoxy [20] specimen by 40 
percent, the corresponding correction factors decrease or increase 1.3 and 2.0 percent, 
respectively. The same insensitivity of the correction factors to shear moduli for 90° and 
0790° specimens is also shown in Fig. 15. 

Simulated Gage Readings 

Sullivan [14] stated that "In a properly loaded Iosipescu specimen, the two strains at 
±45 deg should be equal in magnitude and opposite in sign to a reasonable 
approximation." Morton et al [7] and Lee and Munro [18] showed recently that this 
condition cannot be achieved experimentally for 0° specimens. Lee and Munro tried to 
achieve Sullivan’s condition by adding shims to the specimen and having the fixture 
contact surface carefully ground, but the effort proved to be of no effect. The inequality 
of strains in the ±45° gages can be attributed to the presence of transverse normal strains 
in the gage section and is an inherent property for 0° specimen. The simulated individual 
gage readings for 0°, 90° and 0790° graphite-epoxy specimens are shown in Fig. 16. The 
magnitude of the predicted strains in the simulated ±45° gages are approximately equal 
for 90° and 0790° specimens, see Fig. 16. However, for the 0° specimens, the strain 
recorded by the +45° gage is larger in magnitude than the strain measured by the -45° 
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gage. The absolute value of the ratio of tensile to compressive strain, as determined from 
the ±45° gages, is 0.6 for graphite-epoxy and 0.8 for Kevlar and glass-epoxy specimens. 


Even though the ±45° gages in the 0° specimen do not measure strains of equal 
magnitude, the shear strain calculation is still valid provided strains in the test section are 
constant. The strains measured by the gages can be expressed in terms of strains in the 
specimen coordinate system (x,y), 

SG 1 = 0.5 (y x 4 / + e x 45 +£ , 45 ) (5) 

SG 2 = 0.5 (-Yxy 5 + e x 45 +e y 45 ) (6) 

Y 12 = SG 1 - SG 2 = 0.5 (Y x 4 > 5 + Y xf +e x 45 - £ x 45 + e y 45 - £ y 45 ) (7) 

where SG's represent strain gage readings, as shown in Fig. 2, and Yxy, £x and e y 
represent shear strain, longitudinal and transverse normal strains averaged over the gage 
length with superscripts ±45 indicating the gage directions. Under pure shear, the 
normal strains are zero and the shear strain can be obtained by 

Y 12 = SG 1 - SG 2 = 0.5 (Yxy 5 + Yxy 5 ). (8) 

For 90° and 0790° specimens, the normal strains (e x 45 . e y 45 , Ex' 45 and £y‘ 45 ) in the gage 
regions are negligible and equation (7) reduces to equation (8). However, for the 0° 
specimens, the longitudinal normal strain, e x , is negligible but the transverse normal 
strain, e y , is significant, about 9.5% and 24.4 % of Yavg for Kevlar-epoxy and graphite- 
epoxy, respectively; therefore SGI and SG2 would not measure equal and opposite 
strains. However, if the transverse normal strain distribution is uniform in the region 
where the strain gages are located, the influence of the transverse normal strain field on 
the two gages is the same even when separate strain gage rosettes are used and the 
calculation of the shear strain would reduce to equation (8). That is, if the transverse 
normal strain field is uniform in the gage section, the shear strain calculated from the two 
gage readings in the 0° specimen is a valid shear strain and can be used in calculating the 
shear modulus. Lee and Munro’s suggestion [ 1 8] that the 0° specimen should not be used 
for measurement of the in-plane shear properties does not appear to be correct, if it is 
based solely on the inequality of the ±45° strain gage readings. 
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I nnri-Point Effect 


The W2 fixture distinguishes itself from the antisymmetric four-point-bend method 
by application of uniform displacement over a large section of the specimen edge and 
thus avoids the problem of local crushing resulting from concentrated loads. One 
consequence of using the W2 fixture is the region of load introduction is a function of the 
applied load and specimen ply orientation. Employing the proposed analysis method the 
load distribution varies as a function of applied load. For small applied loads the load 
distribution is near uniform whereas for applied loads where large deformations occur the 
load distribution would approach the concentrated load case. For illustration purposes, a 
uniformly distributed load and a concentrated load case is used to demonstrate their 
effects on shear strain distribution. Fig. 17. Though these two load cases may be extreme 
the results depicted in Figs. 18-20 show that location of the load points has a profound 
effect upon the shear strain distribution along the notch axis for the 0° specimen. For 
graphite-epoxy, the deviation of shear modulus correction factors from the model in Fig. 
6a could be as large as 10%. The deviation may be as large as 7% and 5% for glass- 
epoxy and Kevlar-epoxy respectively. But for the 90° specimen, the load-point effect is 
negligible for graphite-epoxy and is as low as 2% for Kevlar-epoxy and glass-epoxy for 
these two load models. Hence it is possible that the shear modulus obtained 
experimentally from 0° specimens will have large variations. 


The W2 fixture has load points which are farther away from the notch axis as 
compared to those of the W1 fixture, and the induced transverse normal strains in the 
gage section for the 0° specimen are smaller. Hence it might be suggested that if the 
fixture were redesigned such that the load points were even farther away from the notch 
axis, a better shear distribution might be obtained along the notch axis for the 0 
specimen. However, to obtain shear force P in the test section, the equivalent 
concentrated load, as shown in Fig. 17a, is Pa/(a-b) and Pb/(a-b) for the inner and outer 
load pairs, respectively. If the load points are located as shown in Fig. 17b, where b is 
large, the equivalent concentrated forces for the inner force couple are about twice the 
magnitude of those in Fig. 17a and the equivalent concentrated forces for the outer force 
couple are about six times larger in magnitude than those in Fig. 17a. Thus if the load 
points are too far away, the concentrated forces would become so large that local 
crushing of the specimen may occur. It is shown earlier that the transverse normal strain 
in the test section for the 0° specimen would not affect the calculation of shear modulus 
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because the distribution of transverse normal strain is uniform in the gage section. So, 
even a redesigned W2 fixture will not provide pure shear in the test section, it is not 
necessary to move the loading faces of the fixture farther away from notch axis. 


CONCLUSIONS 

Linear finite element analyses of Iosipescu specimen tested in modified Wyoming 
fixture have been performed and evaluated. It is suggested that contact length between 
the specimen and fixture, constraint type and compatibility should be considered in the 
finite element modeling. To better model the fixture-specimen kinematics, an iterative 
analysis procedure was applied. It was found that the finite element modeling was not 
affected by the material systems (such as graphite-, Kevlar- and glass-epoxy) but are 
influenced by the elastic stiffness due to different fiber orientation (0°, 90° and 0 /90 ) 
producing different deformation of the specimen in the fixture. For the 0° specimen, it is 
shown that the pure shear state in the test section was not achievable. However, the shear 
strain calculation is not affected by the uniformity of normal strain distributions. It is 
shown that correction factors are needed for the calculation of shear modulus and that the 
correction factor for the 0° specimen is strongly dependent on the finite element model 
used, but not for 90° specimen. 
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Figure Captions 

Fig. 1. Original Iosipescu (a) and (b) shear fixture and (c) specimen. 

Fig. 2 Modified Iosipescu shear (a) fixture and (b) specimen. 

Fig. 3. Deformation characteristics of Iosipescu shear specimen. 

Fig. 4. Finite element model of Iosipescu specimen. 

Fig. 5. Initial displacement (applied and constrained) distribution on 
Iosipescu specimen. 

Fig. 6. Distributed loads on the specimen sides due to imposition of uniform displacement 
for (a) 0°, (b) 90°, (c) 0790° specimens. 

Fig. 7. Normal and shear strains at x=0, normalized with respect to average shear strain, 

~ for (a) graphite-epoxy, (b) Kevlar-epoxy, (c) glass-epoxy 0° specimen. 

Fig. 8. Normal and shear strains at x=0, normalized with respect to average shear strain, 
for (a) graphite-epoxy, (b) Kevlar-epoxy, (c) glass-epoxy 90 specimen. 

Fig. 9. Normal and shear strains at x=0, normalized with respect to average shear strain, 
for (a) graphite-epoxy, (b) Kevlar-epoxy, (c) glass-epoxy 0790 specimen. 

Fig. 10. Typical shear strain contours for (a) 0°, (b) 90°, (c) 0790° graphite-epoxy 
Iosipescu specimen. 

Fig. 11. Transverse normal strain contours for (a) 0°, (b) 90° graphite-epoxy 
Iosipescu specimen. 

Fig. 12. Longitudinal normal strain contours for (a) 0°, (b) 90° graphite-epoxy 
Iosipescu specimen. 

Fig. 13. Shear strain distribution at x=0 for (a) 0°, (b) 90° graphite-epoxy, 

Kevlar-epoxy and glass-epoxy specimens. 

Fig. 14. Correction factor vs.material orthotropic ratio. 

Fig. 15. Sensitivity of the correction factors to shear moduli for graphite-epoxy Iosipescu 
specimen. 

Fig. 16. Normalized strains of ±45° gages from finite element analysis for (a) 0°, (b) 90°, 
(c) 0790° graphite-epoxy Iosipescu specimen. 

Fig. 17. Force couple finite element models of Iosipescu specimen tested in Wyoming 
" fixture, (a)concentrated load model, (b) distributed load model. 

Fig. 18. Load point effect on the shear strain distribution at x=0 for (a) graphite-epoxy, 
(b) Kevlar-epoxy, (c) glass-epoxy 0° specimen. 

Fig. 19. Load point effect on the shear strain distribution at x=0 for (a) graphite-epoxy, 
(b) Kevlar-epoxy, (c) glass-epoxy 90° specimen. 
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pja. 20. Load point effect on the shear strain distribution at x=0 for (a) graphite-epoxy, 
(b) Kevlar-epoxy, (c) glass-epoxy 0790° specimen. 
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Fig. 1. Original Iosipescu (a) and (b) shear fixture and (c) specimen. 
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Fig. 2. Modified Iosipescu shear (a) fixture, (b) specimen 
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Fig. 3. Deformation Characteristics of Iosipescu shear specimen 



v= 5 



Fig. 5. Initial displacement (applied and constrained) 
distribution on Iosipescu specimen. 
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Fig. 7. Normal and shear strains at x=0, normalized with respect to average shear strain 
for (a) graphite-epoxy, (b) Kevlar-epoxy, (c) glass-epoxy 0° specimen. 
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Fig. 10a Typical shear strain contours for 0° graphite -epoxy Iosipescu specimen. 
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Fig. 10c Typical shear strain contours for 0°/90° graphite-epoxy losipescu specimen. 
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Fig. 1 la Transverse normal strain (e y ) contours for 0° graphite-epoxy Iosipescu specimen. 
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Fig. 12a Longitudinal normal strain (e^) contours for 0° graphite-epoxy losipcscu specimen. 
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Fig. 13. Shear strain distribution at x=0 for (a) 0°, (b) 90°, (c) 0790° graphite-epoxy 
Kevlar-epoxy and glass-epoxy specimens. 
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Fig. 15. Sensitivity of the correction factors to shear moduli for 
graphite-epoxy Iosipescu specimen. 















Pb_ Pa_ 

a-b i a-b 



Fig. 17. Force couple finite element models of Iosipescu specimen 
tested in Wyoming fixture, (a)concentrated load model, (b) 
distributed load model. 


























Fig. 20. Load point effect on the shear strain distribution at x=0 for (a) graphite-epoxy, (b) 
Kevlar-epoxy, (c) glass-epoxy 0°/90° specimen. 






